Stochastic nonlinear differential equation generating 1/f noise 
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Starting from the simple point process model of 1/f noise we derive a stochastic nonlinear 
differential equation for the signal exhibiting 1/f noise in any desirably wide range of frequency. 
A stochastic differential equation (the general Langevin equation with a multiplicative noise) that 
gives 1/f noise is derived for the first time. The solution of the equation exhibits the power-law 
distribution. The process with 1/f noise is demonstrated by the numerical solution of the derived 



equation with the appropriate restriction of the diffusion of the signal in some finite interval 
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The power spectra of a large variety of systems ranging widely from astrophysics and 
technology to sociology and psychology at low frequencies have 1/f behavior, i.e., the power 
density S(f) is inversely proportional to the frequency / . 1/ / noise, 

also known as flicker noise, is intermediate between white noise (no correlation in time, 
S(f) ~ 1//°) and the Brownian motion (no correlation between increments, S(f) ~ 1/f 2 )- 
Simple procedures of integration or differentiation of such fluctuating signals do not yield the 
signal exhibiting 1/f noise. Most of the 1/f noise models are specialized or complicated. 
This makes the problem of omnipresent 1/f noise one of the oldest puzzles in the con- 
temporary physics. In contrast to the Brownian motion generated by the linear stochastic 
equation, simple systems of differential, even linear stochastic equations generating signals 
with 1/f noise are not known. 

The purpose of this paper is the derivation of a nonlinear stochastic differential equation 
(generalized Langevin equation for the signal) generating signal with 1/f noise. The stochas- 
tic differentia l eq uation is obtained from the point process model of 1/f noise, analyzed in 
Refs. 3, II, H H, Q, Ejj, 0. Such a method enables one to obtain various stochastic 
differential equations, starting from different point processes and generating stochastic sig- 
nals with different slopes of the power density. Analysis of the concrete physical models 
and application of the derived nonlinear stochastic equation for modeling of the specific 
observable processes are out of the scope of this paper. 



We start from the point process model recently proposed and analyzed in Refs. 0, 
, 00]. The signal in the model consists of pulses or series of events, 
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I(t) =a^25(t-t k ). 



(1) 



Here 5(t) is the Dirac delta function, {t k } is a set of the occurrence times at which the 
particles or pulses cross the section of observation and a is a contribution to the signal of 
se or particle. The power spectral density of the point process (JTJ) may be expressed 
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where T is the observation time and 

A(k; q) = t k+q - t k 



k+q-1 
l=k 



(3) 




FIG. 1: Power spectral density of the point process, described by the equations © and (JSJ). 
Parameters used are <r = 0.01, T m ; n = 10~ 5 , and 

''"max — 1- The dashed line represents the spectral 
density, calculated according to equation S(f) = 1/ f. 

is the difference between the pulses occurrence times tk+ q and t^. Here the brackets (. . . ) 
denote the averaging over the realizations of the process and = t^+i — is the interevent 



time. In the model 
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14L Il5l Il6| the interevent time of the signal stochastically 



diffuses about some average value and the process has been described by an autoregressive 
iteration with a very small relaxation. Here we will consider the stochastic point process 
described by the recurrent equations 



tk+1 — tk + T k, 
Tk+1 = r k + oe k 



(4) 
(5) 



with the appropriate boundary conditions, restricting the diffusion of in the finite interval 
[r-min , T max ]. In Eq. (jSJ) Sk are normally distributed uncorrelated random variables with a zero 
expectation and unit variance, i.e., a white noise, and a is a standard deviation of the white 
noise. 
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The signal (JH) generated according to Eqs. (jU and ©, depending on the parameter a 
and the interval [r min ,r max ], exhibits 1/ f noise in any desirably wide range of frequency. 



According to the general theory 
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16^ the power spectral density of 



such point process for / < r m£ [ x and r min — > may be estimated as 

m~£-- r (6) 

max J 

The spectrum obtained from the numerical solution of Eqs. (jlj) and (0) with reflective 
boundary conditions at r min and r max is shown in Fig. [TJ We see that the considered point 
process gives 1/f noise in a wide range of frequencies. 

It is the purpose of this paper to derive a stochastic differential equation for the signal, 
the solution of which exhibits 1/f noise. For this purpose we rewrite Eq. © as a differential 
Ito stochastic equation interpreting k as a continuous variable, i.e., 

t = •«*>■ < 7 > 

Here is a Gaussian white noise satisfying the standard condition 

mm')) = s(k-k'). (8) 

Then we rewrite Eq. (JZj) using the occurrence time. Transition from the occurrence number 
k to the actual time t according to the relation dt = r^dk yields equation 

The si gna l averaged over the time interval T/~ according to Eq. (0) is x = ajr^. The 



standard 



17| transformation of the variable from r to x = a/r in Eq. Q results in the 



stochastic differential Ito equation 

+>*««>■ < 10 ' 

Equation (|10|) can be rewritten in the form that does not contain any parameters. Intro- 
ducing the scaled time 

t s = —t (11) 



a 

we obtain from Eq. (fTU|) an equation 



^ = x* + x^i{t s ). (12) 



The steady state solution of the stationary Fokker-Planck equation with the appropriate 
reflective bounda ry c onditions and a zero flow obtained from Eq. (|T2"j) according to the 
standard method [17[ is of the power-law form, 

P(x) = (13) 



where C has to be defined from the normalization. 

The power-law distribution of the signals is the phenomenon observable in a large variety 
of processes from earthquakes to the financial time series 0,0,0,0]- Therefore, our model 
of 1/f noise is complementary to the models based on the superposition of signals with a 
wide-range distribution of the relaxation times resulting into the Gaussian process j^j]. 

Because of the divergence of the power-law distribution and requirement of the station- 
arity of the process, the stochastic equation (J12j) should be analyzed together with the 
appropriate restrictions of the diffusion in some finite interval x m i n ~ £ ~ x max . Such re- 
strictions may be introduced as some additional conditions to the iterative solution of the 
stochastic differential equation. The similar restrictions, however, may be fulfilled by intro- 
ducing some additional terms into Eq. (fT2"j) . corresponding to the restriction of the diffusion 
in some "potential well". According to the general theory the exponentially restricted 
diffusion with the distribution density 

generates the stochastic differential equation 



dt„ 2 V x 



~n—4 



+ x i + x 5/2 £(t s ). (15) 



Here n is some parameter. 

Equations (j!2)l and ()15j) are the main result of this letter. Since the point process (J3J) and 
(jSJ gives the signal with 1/f noise, the signal obtained from Eqs. (fT2"|) and (j!5)l also should 
give 1/f noise in some frequency interval. When x max — > oo, from Eq. © we can estimate 
the power spectral density as 

S(f) - x 2 min j . (16) 
Such a conclusion is confirmed by the numerical solution of Eq. (|15|1. 
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FIG. 2: Distribution density P{x) obtained from the numerical solution of Eq. (|18|) . The dashed 
line represents the distribution density calculated from Eq. dJ). Parameters used are x m \ n = 1, 
a^max = 10 3 , n = I, and n = 0.1. 



We solve Eqs. ()12|) and (|15|) using the method of discretization. When the variable step of 
integration is At s = hi, the differential equation (fTK|) transforms to the difference equation 
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hi Si. 



(17) 



We can solve Eq. (J 17)) numerically with the constant step, hi = const, when = £j + h. 
However, one of the most effective method of solutions of Eq. ()17)) is when the change of the 
variable Xj in one step is proportional to the value of the variable. We take the integration 
steps hi from the equation x^ 2 \fh~i = KXi with k<1 being a small parameter. As a result 
we have the system of the equations 



x i+l 



H+l 



.2 



1 + - 



n x 



X; 



X; 



x. 



(18) 



U + 



K 



r 



The distribution density P(x) of the variable x, obtained using Eq. (fTHj) is shown in Fig. 




FIG. 3: Power spectral density of the signal, obtained from Eq. (|18|). The dashed line represents 
the spectral density, calculated according to equation S(f) = l/(2-/r/). The parameters used are 
the same as in Fig. |2 

We see that our method of solution gives good agreement with the power-law distribution 
(H3J) in the interval x min < x < x max . 

The power spectral density S(f) calculated according to Eqs. (fT%|) is shown in Fig. |31 
Figure El shows that equation (fL3|) indeed gives a signal exhibiting 1/f noise in a wide 
frequency interval. 

In summary, for the first time we derived a stochastic differential equation for the signal 
exhibiting 1/f noise in any desirably wide range of frequency. The distribution density of 
the signal is of the inverse cubic power-law. The numerical analysis of the obtained equation 
shows that the signal indeed exhibits 1/f noise and power-law distribution. 

The research described in this publication was supported in part by the Lithuanian State 
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